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Abstract 

We produce the simple formula for the radius, R, of a cold body: 

where y 3 = M/M p , M is the total mass, J = y/l - y A jy\ h and I = J/ 1 + f ^1 - 

y 3 h = M c h/Mp where M c h is Chandrasekhar's limiting mass. The density of our liquid/solid matter 
at low pressure defines the constant p . For bodies with non relativistic electrons I — J — 1. 

M p is constant and gives the mass of the planet of maximum radius. M c h, M p and p axe all defined 
in terms of fundamental constants, in particular the ratio M p /M c h is basically given by the three halves 
power of the fine structure constant. 

The aim of the discussion is to emphasise physical principles and to connect the masses derived to the 
fundamental constants, so mathematical niceties are sacrificed and replaced by interpolations between 
simple exact limits. 



1 Introduction 

The virial theorem shows that a star deprived of any energy sources will shrink and its kinetic energy will 
grow which normally implies it will get hotter. In the 1920's Eddington (1926) was puzzled as to how a star 
could ever cool down and end its life. This problem, together with the mystery of white dwarfs such as Sirius 
B, was beautifully solved by R. H. Fowler (1926) who realised that the electrons would become degenerate, as 
in a metal, so the star could then be supported by the degeneracy pressure which is an inevitable consequence 
of the zero-point energy combined with the Pauli exclusion principle. 



Even bodies of zero temperature have such a pressure which supports the white dwarfs against gravity. 
Stoncr (1929) made early models of white dwarfs to compare their densities with the theory but Anderson 
(1929) rightly criticised him for not allowing for the relativistic motions of electrons forced to very high 
densities. He showed that when this was allowed for, there would be a limiting mass beyond which no cold 
white dwarf models could exist. His estimate of the limiting mass was not very good, but Stoner (1930) then 
modified his own calculations and calculated the limit for homogeneous models which gave an answer good 
to 10%. Meanwhile, Chandrasekhar (1931) had developed the theory independently starting from Fowler's 
work. He realised that the limiting configuration would be an n = 3 polytrope and gave the first accurate 
determination of the Chandrasekhar limit. 

More detailed models of white dwarfs were given by Salpeter (1967) who derived conditions for their so- 
lidification and determined properties of the lattice. Shapiro and Teukolsky in their book give a definitive 
account of how all this extends to neutron stars. 
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Here our aim is much more modest. To give a rough approximate formula for the radius of any cold body as 
a function of its mass. We illustrate how the balance between degeneracy pressure, the electrical attraction 
of electrons to atomic nuclei, and gravitation lead to a mathematical formula that describes the mass-radius 
relationship for cold bodies. This holds from asteroids, through planets and brown dwarfs to white dwarfs 
and neutron stars. Understanding rather than accuracy is our aim so we shall use interpolation formulae 
where they are simpler, nevertheless we shall ensure accuracy in important limiting cases. Thus we aspire 
to follow the fine tradition set in Salpeter's 1967 lectures. 

At masses below that of Saturn, the degeneracy pressure of the electrons within the atoms (vide the Thomas 
Fermi atom) is balanced by the electrical binding. While the gravity serves to hold the whole body in a 
nearly spherical shpae, it is not strong enough to crush the atoms into significantly smaller volumes. Thus 
the material is nearly incompressible and apart from any chemical fractionation the lighter bodies have the 
same density. However, above the mass of Saturn the gravitational potential energy increases to become 
comparable with the electrical one. Above the mass of Jupiter, the gravitation dominates, the atoms are 
crushed together with the outer electrons eventually booming pressure ionised to make a continuous sea. 
In fact, somewhat above Jupiter's mass, the weight of each layer added so compresses the material within 
that the whole body shrinks in radius. The electron sea is not so dense that the zero point energy drives 
the electrons to speeds close to that of light, so they give Fowler's p = Kp 5 ' 3 as an equation of state as in 
n = 3/2 polytropes. The J pdV energy behaves as M 5 / 3 /i? 2 . Although the —GM 2 /R gravitational energy 
rises with a higher power of the mass, the total energy achieves its minimum at the (equilibrium) radius 
with R oc M~ 4 / 3 . This is the brown dwarf white dwarf sequence. However, as M is further increased the 
increasing density of the bodies (p oc M 2 ) leads to zero point agitation with relativistic speeds. This results 
in the pressure density relationship p = K'p 4 ^ 3 with j pdV oc M 4 ! 3 jR. As this is approached the equilibrium 
radius rapidly becomes very small and very sensitive to mass because both gravitational and internal energies 
scale with the same power of the radius. No force balance is possible beyond Chandrasekhar's limiting mass. 

In neutron stars the electrical energy is unimportant; most of the electrons have combined with protons and 
been transformed into neutrons so that neutron degeneracy supplies the pressure. Neutron star radii are 
smaller than white dwarf radii by the factor ^f- so the densities are (m n /m e ) 3 times greater. 

A secondary aim of this article is to express the important masses that arise in astronomy in terms of 
the fundamental physical constants. Thus the ratio of the Chandrasekhar mass to the mass of the planet 
of maximum radiu is essentially the fine structure constant raised to the power —3/2. Such matters are 
of particular importance to the misguided few who seriously entertain theories in which such fundamental 
constants change. 



2 Kinetic Energy of Degenerate Material 

The kinetic energy of an electron is given by 



where the initial 2 comes from the 2 spin states of an electron. If all levels up to momentum p m are filled, 
the kinetic energy of the electrons per unit volume is 




(1) 



The number of energy levels in a box of side L with momentum between p and p + dp is 



2 x 4np 2 dpL 3 h~ 3 , 




ep dp . 



(2) 



While the integral may be evaluated using sh 
the limiting cases and interpolate. 




that begins to obscure the phyi 



sics so we shall evaluate 
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If p m <C m e c then e = ^p 2 /m e so 



£ = —h 3 p 5 m /m e . 
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If Pm ^> m e c then e ~ pc for most electrons, so 

£ = 2nh~ 3 pf n c 



(3) 



(4) 



Looking at (1), an interpolation in terms of the function \/x 2 + 1 — 1 where x = r/^ 1 ^ suggests itself, where 

2 

i] is to be determined. For large x this function gives r]p m /(m e c) while for small x it gives ±r] 2 . 



So the ratio of the non-relativistic case to the very-relativistic case is iv^Tc- Now the ratio of the expressions 
(3) and (4) is § so we choose r\ = § . The resulting interpolation formula exact at both limits is 



£ = inh^plm.c 2 (V* 2 + 1 - l) 



where 



^Pm_ 

5 TO e C 

If the total mass per electron is p e mn then the density is 

.3 l-3 



where 



p = p e m H ^np m h " = pix" 



3 irp e m H , r , ., 
Pi = (I) 5 / [V(w e c)] 



(5) 
(6) 

(7) 
(8) 



thus pi is the density at which the electrons start to become relativistic due to the uncertainty principle. 
rriH is the mass per electron and h/(m e c) is its Compton wavelength. Rewriting (5) in terms of p 



m e c 2 ( rr, 

£ P=%—?—(\/x 2 + \-l 
rriHUr \ 



Now imagine expanding unit mass of this material 

l 



PdV = Pd = -d (£/p) = -d (£/p) jdx dx . 
Since p oc x 3 we find the pressure density relationship 

^5 

P = Pi- 



where 



4 4 5 

5 7T mtc 



Pl 2? 3 h 3 ■ 

In particular the non-relativistic and extreme relativistic limits give P = Krp r , 



(9) 
(10) 

(11) 
(12) 




K r = 



notice P\j p\ = 



Plpl 


5/3 


Plpl 


4/3 


m e 


c 2 







T = 5/3 non-rel 
r = 4/3 extreme rel. 

T = 5/3 
T = 4/3 



(13) 



(14) 
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Returning to (10) we see that whenever £ oc p r then 

P = (T-1)£ . (15) 

We have of course devised our derivation so that these formulae (8), (13), (14) & (15) are exact cf Shapiro 
k Tcukolsky (1983). 

Three limiting density distributions are important 

1. The uniform density which we shall see occurs. 

2. The non-relativistic T = § polytrope of index § = n = ■ 

3. The extreme relativistic polytrope of index 3, T = §. 

To orient ideas it is useful to consider initially the uniform density case. The kinetic energy is given by (9) 
&(7) 

T = §ttR 3 £ = A G M 



(16) 



^/(CoMV3/R)2 + I _ 

where we have written x = C M 1 ^/R and A a = ±f = 3 -&-, 

^ u I u 16 TYlH^e Pi 

Co = (|7r Pl )- 1/3 . (17) 

Although this uniform case is only realistic for very small x, nevertheless the modifications needed for 
polytropic cases occur only via structure constants of order unity. Thus much understanding can be gained 
by studying the problems with T given by (16) and (17) for all x. To show the modifications needed for the 
polytropic cases we now take the case of a polytrope of index n = j^j. Using (15) 



-I 



47rr 2 £dr = nK r R 3 y 4tt (-J p r d (r/R) = nK r M J C n , 



where 

C « = jjQ{j) WR) , (18) 

is a structure constant for a polytrope of index n which is evaluated in the appendix: C = 1, C 3 / 2 = 
1.7501, C3 = 2.2146. Notice that our uniform density case (16) also gives the M(M/R 3 ) r ~ 1 factor for 
both the non-relativistic F = § and the extreme relativistic r = f cases. Following our earlier interpolation 
procedure using 77 below equation (4), we now look for an interpolation suggested by (16) between the T = | 
and r = § polytropes. This is of the form 



AM 



(M 1 / 3 x 2 



1 - 1 



(19) 



where A and £ are chosen to make the two polytropic formulae exact when (M 1 / 3 /R <C 1 and 3> 1 
respectively. This gives 

C = = 0.7903C o ; A( 2 = C 3/2 A £ - 1.750lA o C o 2 . (20) 

Notice that (19) has precisely the same form as (16) but the values of A and ( differ. 



3 Potential Energy 



The gravitational potential energy is —aGM 2 /R where Ritter's formula tells us that for polytropes a = 
(see appendix). However, for planets and small bodies the electrons are held in electrically rather than 
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gravitationally. We must therefore include in the potential energy some estimate of the electrical potential 
energy. The mean number density of electrons is 



n e = Mj (±irR 3 ii e m H ) 



so we estimate the electrical potential energy as 

_^ ^_ 

-1/3 



r(-^-\ 



fj, e m H 



where the first bracket gives the number of electrons in the whole body and (3 is a dimensionless constant to 
be determined. The gravitational potential energy is proportional to M 2 while the electrical is proportional 
to M 4 / 3 at fixed R. Particular interest centres on the mass at which these two contributions are equal. 
Calling it Mi we find with a = § 

/ g 2 \ 3 / 2 /5/3\ 3 / 2 

M 1 = (5/3/3) 3/2 M - 2 — - m H = /x e " 2 2-308 x 10 30 5 m . (22) 



\GmjjJ n \ 3 ) 



/ 2 \ 3 / 2 

Except for a structure factor of order unity M\ is given by \j§^i J , the ratio of the electrical to gravitational 

forces between two protons, raised to the 3/2 power, times ran- Notice that this mass does not depend on 
Planck's constant. In adapting formula (21) to rocky planets & asteroids we have used p a = 3.65 and fi e = 2. 



4 Energy Minimisation Gives Radius 



The total energy is now given by 



E = AM ( VC 2 M 2 / 3 R- 2 + 1-1 



GM 2 



Mi 
~M 



2/3' 



The equilibrium radius has the minimum E for any R so 



therefore 



A ^ M 2/3 R - 



R_ dE__ { 
M dR ~ ' 

2 a CM 1 / 3 



^2 M 2/ 3i? -2 + I ( R 
equation (25) is readily solved for (M 1 ^ 3 R~ 1 and thence for R 

CM 1 / 3 



G (ikf 2 / 3 + Ml' 3 ) 



R = 



where 



q = f ~A ( M2/3 + M ' /3 



Evidently R = when q = 1 so this gives Chandrasekhar's limiting mass. Since M 2 ' 3 is negligible at this 
limit we have (putting a = f for the n = 3 polytropc r = |) 



(23) 

(24) 
(25) 

(26) 
(27) 



M, 



ch 



aG 



3/2 /n \ 3/ 2 



t. \ 3/2 

he \ run 



Gr 



H 



Me 



(28) 



Which demonstrates that the Chandrasekhar limit is essentially the gravitational fine structure constant to 
the power of — f times the mass of the hydrogen atom. Equation (26) when combined with (27) gives the 
mass-radius relationship 

C 2 A M 1 ' 3 I / M \ 4/3 , x 

R=^7 — Wl-I — I , (29) 



aG MVS + M 2/3 



ch 
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where we have simplified the surd taking account of the fact the M c h 3> Mi . 



Equation (29) shows that there is a maximum radius for a cold planet close to Mi i.e., close to Jupiter's 
mass. Also, for small M, R grows like M 1//3 so the density is constant. This is not surprising; the electrons 
are held in by electricity and we are merely placing together electrically bound neutral objects. As more 
mass is piled on the growth in R slows until close to Mi the weight of the overlying material so crushes what 
is beneath that additional material hardly changes R. Beyond Mi the crushing is so great that additional 
mass only serves to make R decrease like M -1 / 3 . This is the normal white dwarf regime which continues 
until the electrons are so confined that they become relativistic, the radius then shrinks precipitously as the 
Chandrasekhar limit is approached. 

The same behaviour is predicted whether we use the crude homogeneous approximation of Stoner (1930) 
via (16) or the inhomogeneous ones via (19) but the detailed numbers will be different. The fact that for 
M <C Mi the bodies grow at constant density shows that they are not centrally condensed. Thus in that 
regime we should use the homogeneous model. However, for the non-relativistic white dwarfs we should use 
the polytrope of index §, V = § and for the relativistic white dwarfs the polytrope of index 3, T = §. Each 
of these three cases is described by a different value of a which is g- 3 -^- with n = corresponding to the 
homogeneous case. As n varies systematically with mass taking the values 0, § and 3 in the low mass, white 
dwarf, and relativistic cases it is not difficult to devise and interpolation formula for a. However there is 
one further problem of this type. We chose £ and A to fit the white dwarf density profiles but we now find 
that homogeneous profiles are appropriate for M <C Mi. Thus we should interpolate not just a but rather 
(C 2 A/a) which should take the value [£ 2 A /(§)] f° r small M and ( 2 A/a with a = for the white dwarf 

regime. In the whole range the average density increases with mass so ^/(M 2 / 3 + Mi /3 ) decreases as M 
increases (see 29). 

The interpolation 

eA = £A L m1' 3 + mW ^ 

aG 3/5G ^ M V3 + . 8 i 6M 2/3^ K 
where K = 1 + |(1 — J 1 / 2 ) and J = (1 — y i y~^) 1 ^ 2 , gives the correct values for M small and M large 

r 2/3, i 3/2 

and ensures that the density increases. We now define M p = M^70.816 = 1.356Mi. Then writing 
y = (M/Mp) 1 / 3 , our interpolation formula (31) becomes 

( 2 A _ 5( 2 A M 2/3 + M 2 / 3 



aG 3G M* /3 (l + y 2 )K 

and our Mass radius relationship is 



(31) 



(| 7rPo ) 1 /3 i? = M V3_^_ J i (32) 



/ v 1/3 / s 1/3 

with y = and y ch = (^) and I = J/K, 



(!W /3 = /3f(^) 2/3 M^feV /3 : < = 



2 



V a o J m e e 

so p is apart from structural constants the density of hydrogen within the Bohr radius, a a , of the nucleus. 
p depends on Ti through a a . 

(2 \ 3/2 

This is the mass of the cold planet of maximum radius. As stated it is essentially the fine structure constant 
to the § power times the Chandrasekhar mass or alternatively the ratio of electrical to gravitational forces 
between two protons, raised to the three halves power, times the hydrogen atom's mass. 
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5 The Mass Radius Relationship 



Our mass-radius relationship is now defined except for the constant (3 that we inserted in (21)) to allow for 
the extreme crudeness of our estimate of the electrical potential energy. This we shall evaluate by fitting our 
formula to the radius of Saturn. Away from the relativistic regime formula (32) takes the even simpler form 

(I W / ^ = MV3_^ = ^ . (34) 

Now both and y~ 2 are proportional to (3 so for a planet of known M and R we solve for (3 using Saturn 

/ \ 1/3 

(3= \JT Z *) - (/V) = 1-343 - -230 = 1.113 , ( 35 ) 

where in spite of appearances each bracketed term is independent of (3 and p = M/ (|7ri? 3 ) . The same 
calculation for Jupiter gives (3 = 1.161 so (3 — 1.137 is a good compromise differing from each by only 2.1 
per cent. This value of (3 yields p = OAWgm cm~ 3 and M p — 6.24 x 10 30 pm. The density of Hydrogen at 
the relatively low pressure (compared with planetary interiors) of half a Megabar is close to 0.54 gm cmT 3 
(Alavi et al., 1995). At much lower pressures solid Hydrogen and liquid Helium have densities of 0.07 and 
0.12 gm cm~ 3 . For Terrestial planets f3 (and hence p\J Z and Mp^ 3 ) takes higher values appropriate to their 
composition. From (34) we see that the planet of maximum radius has a density of 8p . Thus although M v 
is independent of h nevertheless R max depends on H. 

To summarise our basic result is that the radius R of a cold body of mass M is given by 




where the surd and the curly bracketed expression reduce to 1 for non-relativistic white dwarfs and small 
bodies. 

According to (36) the cold Planet of maximum radius has a mass of 3.3 times Jupiter's and its radius 
is 8.65 x 10 9 cms as opposed to Jupiter's radius of 7.0 x 10 9 (after correction at constant volume for the 
eccentricity). Equation (36) gives 7.0 x 10 9 for the radius of a body of mass Mj = 1.899 x 10 30 gm 



6 Neutron Stars 



In the above we have not aimed to treat the mass radius relationship for neutron stars but with suitable 
modification the same principles apply to them. 



I It is the degenerate neutrons that provide the support so m e must be replaced by m n wherever it 
occurs. The electric term is small and irrelevant. The mass per neutron is now m n so p e is replaced 

by 1. In particular under (14) pi/pi — j^c 2 and in (8) pi is bigger by the factor ■ 

II The Chandrasekhar mass is almost unchanged except that we are now interested in p e = 1 rather than 
the normal 2. 

III General Relativistic corrections are no longer so small, so the effective gravity is somewhat stronger. 

IV The forces between neutrons are such that the free neutron gas is no longer a very good approximation. 
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Radius (km) 

Figure 1: The dotted line gives the non-relativistic approximation that omits the surd and the bracketed 
expression. The full line gives the full expression for p Q — 3.65gmcm~ 3 , [i e — 2 appropriate for rocky 
planets and white dwarfs. The lighter line gives the expression for a 25% Helium 75% Hydrogen body or 
a pristine white dwarf. The squares are solar system bodies supplemented by seven white dwarfs of known 
radius and one neutron star. The neutron star line is parallel to the dotted lines for non-relativistic white 
dwarfs but reduced in radius from the pure hydrogen white dwarf by the mass ratio of the neutron and 
electron. 
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If we were to ignore all the refinements in III and IV the same formulae hold provided that p a is replaced 

by (^f-^j Po whenever we consider neutron stars. Thus the radii are smaller, by the factor than the 

white dwarf of the same mass with fi e = 1. The Chandrasekhar limit is almost unchanged at 5.82M . [] In 
practice effects III and IV reduce this to between two and 3.6 solar masses, see Shapiro & Teukolsky. These 
basic points were understood by Baade & Zwicky (1934) writing soon after the neutron was discovered. 
Later they took spectra of the Crab Pulsar without realising that it pulsed thirty times a second, but Baade 
could not interpret the continuous spectrum as, in spite of Schott's (1912) work, Synchotron radiation was 
not known in astronomy at that time. 
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A Ritter's formula and the evaluation of the C r , 



1, $ = o at £ = 0. Let the edge be at £ = £ 



Define the "mass" within £ by /i(£) = £ 6»"£ 2 d£. 
Then _£ 2 |= M andf = £ 2 0", 
also £ 3 ^ = (n + l)£ 2 0«£2|/£ = _(„ + i)| ££. 
Thus 



6 f ]f)n+l pCi nil 

£ 3 ^— d£ = -3/ ^+ 1 £ 2 d£ = -(n+l) / £d M , (Al) 
o <K Jo Jo £ 



1 Using the values of the Fundamental Constants & the Sun's Mass then available Chandrasekhar (1938) got 5.75Af0 
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Now write ijj = 9 + ipi and take ip to obey 



1 d 



e delude 



£>£i 



Furthermore choose so that ip — > at oo. 



Hence 



where 



f°° M 

^i = mi/6 v> = 



/i = m for £ > £i , 

Jo 6 £ 2 0" +1 ^ = jf 1 ^« = Jo 1 W> - = Jo 6 « 

using (2) 

= 2 JT 1^-1 = ^/^ by (1) 



Hence 



which gives Ritter's formula 



5 — n f jid^i /if 



-V 



3 GM ; 
5-rc ~~ R~ 



and /% 2 r +1 ^ = — £ 
Jo 5-n 6 



Now 



W — _ . . 1 — 



i "+ 1 M l 



Mi 



— _(t2 M\ 

v 4 



From tables of polytropes 



n = 3/2 £i = 3.65375 m = 2.71406 
n = 3 Ci = 6.89685 u x = 2.01824 



Hence 

C 3/2 = 1-7501 
C 3 = 2.2146 
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